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Ballistic particles interacting with irregular surfaces are representative of many physical problems 
in the Knudsen diffusion regime. In this paper, the collisions of ballistic particles interacting with 
an irregular surface modeled by a quadratic Koch curve, are studied numerically. The q moments of 
the source spatial distribution of collision numbers /j,(x) are characterized by a sequence of "collision 
exponent" r(q). The measure fi(x) is found to be multifractal even when a random micro-roughness 
(or random re-emission) of the surface exists. The dimensions f{a), obtained by a Legendre trans- 
formation from r(q), consist of two parabolas corresponding to a trinomial multifractal. This is 
demonstrated for a particular case by obtaining an exact f(a) for a multiplicative trinomial mass 
distribution. The trinomial nature of the multifractality is related to the type of surface macro- 
irregularity considered here and is independent of the micro-roughness of the surface which however 
influence the values of a m i„ and a ma x- The information dimension Dx increases significantly with 
the micro-roughness of the surface. Interestingly, in contrast with this point of view, the surface 
seems to work uniformly. This correspond to an absence of screening effects in Knudsen diffusion. 
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I. INTRODUCTION 



\Q • Multifractals appear in a wide range of situations like energy dissipation in turbulent flows , electronic eigenstates 
at metal insulator transition 2j, fluctuations in finance dynamics of human heartbeat Q and many others. 
,— 1 i Recently, multifractal behavior of the first passage time has been studied in the transport of diffusing particles on a 
2D Sierpinski gasket with absorbing and reflecting barriers p| . Our interest is to investigate the multifractal properties 
associated with the collision spectrum of ballistic trajectories of point particles interacting with irregular surfaces. 

The emergence of the concept of fractal geometry has provided an efficient tool to model the influence of strong 
geometrical irregularity on various physical and chemical processes 0. Ballistic trajectories in systems with irregular 
geometry are representative of several physical situations. For example, they correspond to the path of light rays in 
irregular structures when diffraction phenomena can be neglected. They also represent the path of atoms or molecules 
I i in confined vessels at sufficiently low pressure, i.e.; in the Knudsen diffusion regime 0. The interaction of atoms or 
molecules with irregular surfaces plays a major role in heterogeneous catalysis |8( where porous catalysts have often 
very irregular surfaces down to the molecular scale. Due to the geometrical irregularity, the small scale structure may 
confine reactants and increase the interaction with the surface. Xenon nuclear magnetic resonance is used to probe 
porous structures through the interaction of rare gas atoms with the solid surfaces 0. The properties of ballistic 
trajectories in 2D pre-fractal surfaces has already been studied in connection with the chaotic aspects, possibility of 
"open" black body, Knudsen diffusivity, mean interaction, and catalytic efficiency by Santra et al[10j . Coppens and co- 
workers have also investigated the influence of the 3D fractal surface morphology on Knudsen diffusion in nanoporous 
media Multifractality has already been observed in the molecular diffusion on random fractal structures with 
time-dependent random potential |l2j| . 

The interest here is to explore the multifractal aspects of the number of collisions made by 2D ballistic trajectories 
within irregular surfaces in the Knudsen diffusion regime. The spatial distribution of the collision numbers of these 
ballistic trajectories over the slit source is found to be multifractal. 



II. THE MODEL 



The macro-irregularity of the surface is modeled by a quadratic Koch curve of fractal dimension In 5/ In 3. The 
pre-fractal generator is shown in FigQ] (a). Particles are launched from a slit of unit length, placed at an angle <j> 
with the zeroth order of the pre-fractal. The particle trajectories due to specular collisions within the third (v = 3) 
generation of the pre-fractal are shown in Fig|I](6). Solid lines with an arrow at the end represent the trajectories. 
The slit source is represented by a thick line extended from x = to x = 1 placed at <\> = 45° with the zeroth stage 
of the pre-fractal surface, a horizontal line. 
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A multifractal measure is in general related to the distribution of a physical quantity on a geometrical support 
[l3j. Here, the geometrical support is the slit source from which the trajectories are launched and the multifractal 
measure n(x) is the distribution of the relative collision numbers, ratio of the number of collisions made by a particle 
to the total number of collisions made by all the particles, over the slit. The particles leave the slit normally. Each 
trajectory is followed as long as they are within the pre-fractal surface. A particle trajectory is described by the 
successive collisions with the surface. The collision points are determined by solving the linear equations describing 
the particle trajectory and the elements of the pre-fractal surface. The trajectories in Fig^fr) are generated using 
this algorithm. 

For the pre-fractal i^th generation, the slit is divided into 3 y equal segments of length l/i v each. A particle is 
placed at the center of each small part. There are then 3 y particles launched from the slit. The number of collisions 
n(x) made by a particle initiated from the position x is counted and the relative collision number fj,(x), ratio of n(x) 
to the total number of collisions n(x) made by all the particles, is then calculated. The value fi(x) is assigned to 
the slit segment positioned at x. 

Generally, particle reflections with surfaces at finite temperatures are non specular and the angle of reflection, 
or re-emission 9 r is different from the angle of incidence 0j. Due to the interaction with the surface, the angle of 
reflection has a distribution around the angle of incidence. Also the reflection angle could be arbitrary due to the 
micro-roughness of the surface. Non-specular reflection for atoms or molecules may come from the complex nature of 
the particle-wall interaction 14]. At the time of hitting the surface, the particle may stay adsorbed for some time and 
then may desorb, exchanging momentum with the surface. For this reason, the particle trajectories as a whole have 
a random character. In order to mimic the random character of each collision, the angle of reflection 9 r is supposed 
to be related to the angle of incidence 9i as 

9 r = 9i + 59 (1) 

where 69 is a random angle distributed uniformly over a range of angle ±r°. This range describes qualitatively the 
random character of the collisions and T can also be called the surface micro-roughness. In the following, the T values 
are expressed in units of degree. A few non-specular trajectories computed using the above rule, are shown in Fig^c). 

In the case of specular reflections, the reflection angle is not determined from numerical computation of trigonometric 
functions since the elements of the pre-fractal geometry are just parallel to the x and y axes. Consequently, there 
is no numerical error due to truncation of angle computations. The only truncation errors come from the solution 
of the linear equations. Because of this error, there could exist trajectories, besides real ray splitting between close 
trajectories hitting salient corners, which are modified by a "spurious" splitting due to number truncation. To avoid 
spurious beam splitting, the cumulative error on the point of collision due to the numerical precision of the computer 
is calculated. For specular reflections, the initial positions of the particles are chosen in such a way that this error is 
always less than the distance between the point of collision and the nearest corner. For non-specular reflections, those 
trajectories are rejected and the particles are resent from the same position with a new random number sequence. 

It should be emphasized that for <fi = 45°, the number of specular reflections for the particular Koch curve considered 
here can be determined from a simple recursion relation without solving the linear equations for the trajectory and 
the surface elements. The number of collisions up to the third generations of the pre-fractal surface are listed in Table 
1 . In the first generation, all the particles send from the central part of the slit have 3 collisions whereas the particles 
from the other two parts have only 1 collision each. The collision numbers are then in the form of (1,3,1) in the first 
generation. In the next pre-fractal generation, each part of the the slit is again subdivided into another 3 parts. It can 
be seen that for v — 2, the collision numbers (3,9,3) for the particles launched from the central part of the slit is just 
3 x (1,3,1) where 3 is the number of collisions corresponding to the same part of the slit for the previous generation. 
The other two sequences (1,3,1) for the two other parts of the slit are given by 1 x (1,3,1). The trajectories from 
these two parts had a single collision each at the previous generation. The multiplication of the collision numbers 
from the first to the second generation of the pre-fractal is shown in Fig[21 Thus, the number of collisions at any two 
consecutive pre-fractal generations are given as: 



v : ni, n 2 , n 3 , n n . . 

U+l: ni,3ni,ni, n 2 ,3n 2 ,n 2 , n 3 ,3n 3 ,n 3 , ■ ■ •, n n ,3n n ,n n 

A recursion relation for the collision numbers for v + 1 generation then can be written as 

(ni, 712,713)^+1 = n v x (1,3,1) (3) 

where n„ is the number of collisions of the trajectory from the same part of the slit at the vih generation. Using 
the above recursion relation, the exact specular collision numbers for 2> v trajectories at any generation v can then be 
obtained for cf> = 45°. 
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III. RESULTS AND DISCUSSIONS 



The collision numbers and their spatial distribution along the slit are computed for the 7th and 8th pre-fractal 
generation (y = 7 and 8). A number of N tot — 3^ equally spaced particles are launched from the slit. The spatial 
distribution of the collision numbers is studied under two different conditions. Firstly, the distribution of the specular 
collision numbers at 45° angle of incidence for which the exact collision numbers are available, is investigated. Secondly, 
the effect of the surface micro-roughness (T) on the collision spectrum is studied for the same slit angle. In this case, 
the collision spectrums are averaged over 160 samples. The general quality of the numerical results is verified with 
the exact collision spectrum for <fi = 45° and T = 0° obtained from the recursion relation (Eq. |3J). 

Before discussing the multifractality associated with the spatial distribution of the collision numbers, it is important 
to consider their statistical distributions. The probability to have a trajectory with N collisions is measured by 

where njv is the number of trajectories having N collisions and N to t is the total number of trajectories. The distribution 
of Pjv for different T in Fig. [3] It can be observed that most of the trajectories exhibit between 3 and 9 collisions. A 
large number of trajectories then leave the pre-fractal surface making only a few collisions and only a few trajectories 
make a large number of collisions. It is important to note that the roughness T has little effect on the collision 
number distribution except the reduction in the maximum number of collisions with increasing T. The total number 
of collisions is also very weakly dependent on the randomness T. These features have already been observed in the case 
of interaction of ballistic trajectories with pre-fractal pores [Tol ] in the Knudsen diffusion regime where the collision 
number distribution was found to be of the Levy type. However, up to the 8th pre-fractal generation, no simple power 
law distribution of the collision numbers was found. 

To study the multifractal aspect of the trajectory distribution, a collision spectrum is obtained. The collision 
spectrum is a distribution of the relative collision numbers of the trajectories over the slit. The relative collision 
number n(x) of a particle trajectory originating from a position x on the slit and making n[x) collisions with the 
surface is given by 

where the sum is over all the particles. This normalized collision number is the measure under study here. In Fig. 
n{x) is plotted against the initial particle position keeping T = (specular reflection) for v = 7. The value of n(x) 
is represented by impulses at the initial positions x. The spectrum is symmetrical around the central position, as 
expected. Most interesting to note, this spectrum corresponds to a uniform exploration of the irregular surface as 
there is a single collision per segment (see Fig. |2J). The maximum number of collisions made by the central trajectory 
is 3 7 and the total number of collisions made by all the particles is 5 7 (the first generation has 5 collisions) at v = 7. 
Thus, [iraax is given by ^ max = 3 7 /5 7 w 0.03, as found numerically. 

The collision spectrums obtained in the presence of surface micro- roughness, r = 2° (a), Y = 6° (b), and T = 10° 
(c) are shown in Fig. [SJfor the v = 7th pre-fractal generation. Due to micro-roughness, the spectrums become more 
uniform and the value of \x m ax is decreased with respect to the specular value. The value of [i m ax is decreasing with 
increasing T as it is also seen in the collision number distribution (Fig. |3J) that the maximum number of collisions is 
decreasing with T. 

To show the multifractal nature of the distribution n(x), it is necessary to study the scaling of the g-moments of 
the measure over different length scales on the slit. If the measure fJ-{x) is multifractal and the slit is divided into n e 
boxes of size e, then the weighted number of boxes N(q,e) is given by 

iV( g ,e) = £/4«e-^ (6) 

i=l 

where fii is the sum of the relative collision number of trajectories initiated in the ith box. Here r(q) is called the 
"collision exponent". In Fig. the values of r(g), obtained by the box counting method, are plotted against q for 
v = 7 (circle) and 8 (square) with T = (specular reflections). The values of r(q) do not depend on v, the macro- 
irregularity of the surface. Note that, r(0) is « 1 and r(l) is « here. r(0) corresponds to the dimension of the slit 
which is 1 here and r(l) is zero because J2 X = 1- The fact that a sequence of exponents are obtained irrespective 
of the generation v confirms the multifractal nature of the spectrum. Note that the multifractality obtained for the 
specular collision spectrum at 45° angle of incidence corresponds to a uniform exploration of the irregular surface 
(Fig. |2l • Consequently, a measure defined over the irregular surface instead of the slit will not exhibit multifractality. 
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A check of the validity of our numerical computations can be made for the case <f> = 45° with specular reflections. 
In that case there exists a trinomial multiplicative cascade |15). A comparison between the collision spectrum and 
the multiplicative cascade is made here. At the zeroth generation of the cascade k — 0, a uniform distribution of unit 
mass over an interval [0,1] is considered. In the next generation k = 1, the unit mass is distributed over three equal 
intervals [0, 1/3], [1/3,2/3], and [2/3, 1] with weighted probabilities p Q — 1/5, pi — 3/5, and p 2 — 1/5 similar to the 
collision number distribution over the slit at ip = 45°. At k = 7 generation, the mass distribution will be exactly the 
same as the collision number distribution over the slit obtained for v = 7 generation of the pre-fractal surface. The 
measure over the unit interval remains always conserved in the iteration process since ^ pi = 1 . The mass exponent 
T m (<l) can be calculated in terms of pi as 

k^oo in 3 h In 3 

The solid line in Fig. H3 represents the mass exponent r m (q) for this trinomial mass distribution. With no surprise, 
there exists an excellent agreement between the exponents for the specular collision number distribution at <f> = 45° 
and the trinomial mass distribution. 

The effect of micro-roughness (T) on the collision exponents is discussed now. In Fig. [7| r(q) is plotted against q 
for r = 2° (circle), 6° (square), and 10° (triangle). The solid line represents the collision exponents for T = 0°, the 
specular reflections. At large positive q, the micro-roughness T has a significant effect on r(q) with respect to V = 0°. 
It is expected because the spectra change with T. Again a sequence of collision exponents r(q) is obtained whatever 
the micro-roughness T. This means that the spectrum are still of multifractal nature. 

To determine the nature of the multifractality, the dimensions/ (a) are obtained as usual as a function of the 
Lipschitz-Holder exponent a through a Legendre transformation 



.15] 



a(q) = -—T(q), f(a)=qa(q)+r(q) (8) 

aq 

of the sequence r(g). The curve f(a) for specular reflexions is discussed first. In Fig. [8] it is represented by circles. 
The solid line represents f(a) for the trinomial mass distribution. Naturally, there is a good agreement between f(a) 
of specular collision distribution and the trinomial mass distribution. Since the slit here is a one dimensional object, 
the f(a) curves are always < 1 and f max (ce) is equal to 1. One observes two parabolas in the f(a) curve. This means 
that the multifractal is of trinomial nature. The values of a m i n and a max for specular reflections at <f> = 45° can 
be predicted. The increment in the measure over a length £ to £ + 5 is = /i(£ + S) — /i(£). The Lipschitz-Holder 
exponent a(£) is defined as /i£ = S a ^'. After nth iteration, 8 = 1/3™ and o/(£) is given by 

(f\ - _ | Inpp + (!-£) ln(l - p - p 2 ) + | lnp 2 , . 

a( ^>~ \nS ~ In 3 [ ' 

where po — P2 = 1/5 and p\ = 1 — po — p 2 = 3/5 for the case considered here. a m i n and a max correspond to £ = and 
£ = 1 respectively and they are obtained as a m in = — ln(3/5) / In 3 « 0.465 and a max = — ln(l/5) / In 3 m 1.465. These 
two extreme values are marked by two crosses on the a-axis. It can be seen that a m i n and a max obtained numerically 
for the specular collision spectrum at <j> = 45° are very close to the predicted values. Note that, as expected, the 
whole f{a) versus a is independent of the generation v of the pre-fractal surface. This is shown in Fig. |SJ where the 
data for v = 7 (circles) and v = 8 (squares) collapse onto the same curve. It is also interesting to notice that a max 
is exactly equal to the dimension In 5/ In 3 of the pre-fractal and a m i n — a max — 1. These features are specific of the 
prefractal geometry and the angle of incidence considered here. 

Always in the case of specular reflections, but for other angles of incidence like cf> = 30° and 60°, it is found that 
the f(a) versus a curves also consist of two parabolas. However, in both cases, the values of a m in and a max keep on 
changing with the generation v upto 8th generation. This is attributed to the fact that the multiplicative process in 
these cases is more involved than in the case of 45° incidence. There should exist for these cases a cascade of different 
sequential values of po, Pi and p 2 so that the asymptotic multifractal measure is still not reached for v — 8. 

The effect of micro-roughness of the surface on the /(a) curve is discussed now. In Fig. |3 f{a) versus a curves 
are plotted for r = 2° (circles), 6° (squares), and 10° (diamonds) at <f> = 45° and v = 7. The solid line represents the 
f(a) curve for specular collisions. The triangles represent the results for an hypothetic "totally" diffusing surface for 
which re-emission would be isotropic, i.e.; occuring randomly in any direction between —90° to 90° with respect to 
the normal, independent of the angle of incidence. This is equivalent to the T — > 90° limit. There arc two important 
things to note. First, the micro-roughness (or random character of the particle re-emission) of the surface is unable 
to destroy either the multifractality or the trinomial nature of the multifractality even at the maximum possible 
value of r, corresponding to the "totally diffusing surface" . Second, the values of a m i n depends strongly on the 
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degree of surface roughness T whereas a max remains almost constant. The value of a m i n increases as a function 
of r. This is expected. As T increases, the value of \i m ax decreases, the spectrum becoming more uniform as seen 
in Fig. and the smaller value of [i max corresponds to the larger value of a m ;„. The value of a max could also be 
estimated from the knowledge of the minimum measure [i m i n . For a single realization, the total number of collisions 
is ks 10 5 . Thus, fimin ~ 1/10 5 (for a trajectory having a single collision) and the value of a max is then given by 
In Umin I In 5 — ln(l/10 5 )/ ln(l/3 7 ) « 1.497 very close to what is obtained here. Since the value of fimin is almost 
independent of the surface roughness (total number of collisions almost independent of T) in this collision process, 
the value of a max remains unchanged. 

From these results, it is found that the multifractality as well as the trinomial nature of the multifractality asso- 
ciated with this collision process remains unaltered irrespective of the surface micro-roughness. This is because, the 
multifractality of the collision spectrum is the outcome of a multiplicative cascade of the spatial sub-distributions 
which are independent of the various situations considered here. The trinomial nature of the multifractal is due to the 
type of pre-fractal chosen here to model the surface macro-irregularity: a quadratic Koch curve, each element being 
divided into three elements in the next generation. A different pre-fractal generator, i.e.; a different model of surface 
macro-irregularity will then lead to a different multifractal behavior. 

The information dimension |l6j| is the value of f(a) when f(a) = a. The information dimension Dj is found to 
be Di ss 0.854 for 45° angle of incidence with T = and at v = 7. This indicates that the measure is supported 
by a minor fraction of the particle trajectories. This also illustrates the fact that most of the trajectories interacts 
weakly with the fractal surface. Interestingly, the information dimension increases with the random character of the 
particle-surface interaction. Indeed for T = 2°, 6° and 10° the values of information dimensions are Dj « 0.877, 0.934 
and 0.952 respectively. For the diffused surface Dj w 0.964. This corresponds to more and more uniform collision 
number distribution on the slit. 

It is interesting to comment on the uniform exploration of the surface, each element receiving the same number 
of particles at least for an angle of incidence of 45 degrees. This is analogous to the results obtained by Andrade et 
al [TtJ through a molecular dynamics study of Knudsen diffusion in the same geometry. In this last case, it is the 
first collision that is uniformly distributed over the perfectly absorbing surface whereas here, the surface is perfectly 
reflecting. The uniform exploration found in these both cases constitutes a qualitative evidence of the absence of 
screening in Knudsen diffusion. It can be considered as a remarkable fact that, seen from the source of particles, the 
collision numbers are strongly non-uniform, while seen from the surface the exploration is uniform. This may have 
interesting consequences for catalyst deactivation which should also be uniformly distributed. 

IV. CONCLUSION 

The collision number distribution and their spatial behaviour are obtained for 2D ballistic particles interacting 
with a pre-fractal curve of dimension In 5/ In 3. The trajectories considered here correspond to the Knudsen diffusion 
regime. The micro-roughness of the surface, or the random character of the particle-surface interaction, is included 
through a randomness parameter T in the reflection angle. It is found that most of the particles make only a few 
collisions whereas a few trajectories exhibit a large number of collisions. A multifractal measure of the spectrum 
is defined in terms of the relative collision number over a linear slit from which the trajectories are launched. A 
sequence of "collision exponents" r(q), is determined by the box counting method and it is found that the sequence 
depends on the moment q. Thus, the spatial distribution of the collision numbers is a multifractal spectrum. The 
fractal dimensions f(a) depend on the random character of the particle re-emission (r). Irrespective of the values 
of r, the plot of /(a) versus Lipschitz-Holder exponent a consist of two parabolas. The spectrum is then that of a 
multiplicative trinomial multifractal independent of surface micro-roughness. The information dimension Di of the 
measure is found to be smaller than 1 and increases with the random character of the particle surface interaction. 
Different type of surface macro-irregularity would lead to different multifractal behavior. Notably, the surface however 
works uniformly, indicating an absence of screening in the Knusdsen diffusion regime. 
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TABLE I: Number of specular collisions made by ballistic trajectories incident at 45° and interacting with a quadratic Koch 
curve of dimension In 5/ In 3 at different generations v. The number of collisions at a generation v + 1 is given by n v x (1, 3, 1) 
where n v is the number of collisions in the previous generation for the same trajectory. 



(a) Pre-fractal generator 




(b) (c) 



FIG. 1: Schematic representation of the collision process at the third generation of the pre-fractal surface, a quadratic Koch 
curve of dimension In 5/ In 3. The generator of the pre-fractal is shown in (a). The thick lines, extended from x = to x = 1, 
in (b) and (c) represent the slit. The slit angle <j> is 45° here. In (b), the trajectories of particles with specular reflections 
form three different position on the slit are shown. The central trajectory makes 3 3 = 27 collisions and the other two make 9 
collisions each. The trajectories of non specular reflection are shown in (c). These trajectories are generated for r = 2, 6 and 
10 for the same initial positions x and slit angle <f> of (b). 




(a) (b) 

FIG. 2: Multiplication of specular collision numbers from v = 1 generation of the per-fractal (a) to v = 2 generation of the 
per- fractal (b) at <j> = 45° . Note that the collision numbers are uniformly distributed over the pre-fractal surface, each segment 
has one collision. 
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FIG. 3: Collision number distribution Pn versus N. The distributions are given for T = (circles) specular reflection, 2° 
(squares), 6° (triangles), and 10° (diamonds) at 45° slit angle. The distribution shows that most of the particles make few 
collisions while a few trajectories exhibit a very large number of collisions. Note the existence of a tail in the distribution. 
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FIG. 4: Plot of the relative collision numbers fi(x) versus x, initial position on the slit, for F — 0° specular reflections. This 
"specular" spectrum is symmetric around the central position as it is expected. 
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FIG. 5: Plot of relative collision numbers fi(x) as a function of the initial position (x) on the slit for (a) T = 2°, (b) T — 6°, 
and (c) r = 10°. Note the general decrease of the values and the increased spreading of the distribution as F increases. 




FIG. 6: Plot of the numerical r(q) against q for different generations v = 7 (circles) and v = 8 (squares) with F = 0. The 
collision exponents remains almost the same. The solid line represents the theoretical r m (q) for the trinomial mass distribution. 
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FIG. 7: Plot of r(g) against q for T = 2 (circles), T = 6 (squares), and r = 10 (triangles) at 45° slit angle. The solid line 
represents r(q) for the specular collisions. The randomness factor T has a significant effect on r(q) for large positive q values. 




FIG. 8: Plot of /(a) against the Lipschitz-Holder exponent a for specular collisions at v = 7 (circles) and v — 8 (squares). The 
crosses represent the values of a m in and a ma i and the solid line represents f(a) for trinomial mass distribution. The curves 
are identical. 




FIG. 9: Plot of /(a) against a for different F at cj> = 45°: circles for r = 2°, squares for T = 6° and diamonds for F = 10°. 
Triangles correspond to the totally diffusing surface and the solid line corresponds to /(a) of the trinomial mass distribution. 



